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Abstract 

We present in this article all Hamiltonian systems in E{2) that are sep- 
arable in cartesian coordinates and that admit a third-order integral, both 
in quantum and in classical mechanics. Many of these superintegrable sys- 
tems are new, and it is seen that there exists a relation between quantum 
superintegrable potentials, invariant solutions of the Korteweg-De Vries 
equation and the Painleve transcendents. 
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I Introduction 



In classical mechanics, an n-dimensional Hamiltonian system is called Liouville 
integrable if it allows n functionally independent integrals of motion in involu- 
tion (including the Hamiltonian), that is 

{H,X,}^0, 

{x,,x,}^o,yij. 

The Hamiltonian H = H(xi, ...x„,pi, ...,pn) and the integrals of motion 
Xi = Xi{xi, ...Xn,Pi, ■■■,Pn) must be well defined functions on phase space 
llUj'). The system is superintegrable if it allows more than n functionally inde- 
pendent integrals, n of them in involution. It is called maximally superintegrable 
if it allows 2n — 1 integrals of motion. The best known superintegrable systems 
in n dimensions are the harmonic oscillator V — ujr'^ and the Coulomb po- 
tential V = J, and they are indeed maximally superintegrable. This may be 
closely related to Bertrand's theorem (Ql O) which states that these are the 
only rotationally invariant systems for which all finite trajectories are closed. 

In quantum mechanics, a Hamiltonian system is said to be integrable if there 
exists a set {Xi\ of n well defined, algebraically independent operators (includ- 
ing the Hamiltonian) that commute pairwise. It is superintegrable if it possesses 
further independent operators, {Yj} that commute with the Hamiltonian. The 
Yj do not necessarily commute with each other, nor with the Xi. 

The independence of operators in quantum mechanics remains to be defined 
rigorously ^1 El El I2HI' Since we are dealing here only with polynomial 
differential operators, we can proceed by analogy with the classical case, keeping 
in mind that a rigorous definition will be needed as soon as we will want to make 
some more general statements. This choice of a definition will be used only 
for discussion purposes, since we will find all potentials that admit third-order 
integrals and all their integrals. The results obtained will therefore hold for any 
definition of the independence of operators. 

Integrable and superintegrable systems, both in quantum and in classical 
mechanics, attracted considerable interest in the last years. Extensive litera- 
ture exists about systems with second-order integrals of motion, either in eu- 
clidian space [SI [3 |S1 El 1^ UHli or in spaces with nonzero constant El 1^ 
or nonconstant curvature El- long as there was no magnetic field in the 
Hamiltonian, the quantum and classical integrals of motion obeyed the same 
determining equations, and therefore quantum and classical integrability were 
very similar. Both properties were related to separation of the Hamilton-Jacobi 
or Schrodinger equations, and also to exact solvability El and generalized sym- 
metries |25| . 

Systems with higher-order integrals have been studied and classified as early 
as 1935 in a well-known paper by Drach. This paper considered classical Hamil- 
tonians in complex Euclidian space. Efforts were made recently to understand 
and classify more completely systems with higher-order integrals in classical 
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[121 177| and quantum mechanics jj^l ■ In spite of these efforts still relatively few 
such systems are known. This paper is the logical sequel of a systematic search 
for superintegrable systems with higher order integrals started in |12j . Here we 
consider two-dimensional real Euclidian space with a one-particle Hamiltonian; 

We request the existence of two additional integrals of motion, one of second 
order in the momenta and the other of third order. 

The condition of existence for second-order integral implies, both in classical 
and quantum mechanics, that the Hamiltonian be separable in cartesian, polar, 
parabolic or elliptic coordinates. In this paper we consider potentials that are 
separable in cartesian coordinates; 

We found all such systems which admit third-order integrals. Quantum 
and classical mechanics will be treated simultaneously, for the conditions of 
existence of integrals of motion, even though not equivalent, are quite similar 
in both cases. 



II Existence of a third-order integral 

In quantum and classical mechanics, the general third-order commuting opera- 
tor 

3 

i+J=0 

can be reduced to a much simpler form, 



^= ^i3k{L3,pipl} + {gi{x,y),p^} + {g2{x,y),py} 

L3 ^xpy - ypx- 

where the Aijk are real constants, and the gt real functions. The bracket is 
the anticommutator. It is not needed in classical mechanics, but in quantum 
mechanics it allows us to get rid of terms with even powers of the pi and to make 
sure the operator is self-adjoint. Furthermore, its use allows us to see clearly 
the relations between the quantum and the classical case. Indeed, it was found 
in | 12j that the requirement that the operator commutes (or Poisson-commutes) 
with the Hamiltonian 

H=^{pl+pl)+V{x,y). 
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implies equations that behave well in the classical limit. Namely, commutativity 
implies 



— giVx + g2Vy — ( flVxxX + flVxxy + faVxyy + fiVyyy 

4 V (2) 
+ 8A3oo{xVy ~yVx) + 2{A2wVx + A2oiVy)j, 

{9i)x = ih{v)Vx + f2{x,y)Vy, (3) 
{92)y = h{x,y)Vx + i.f^{x)Vy, (4) 
{9i)y + (.92). = 2{.f2{x,y)Vx+h{x,y)Vy), (5) 

in quantum mechanics, where 



= -^300Z/^ + ^210Z/^ - ^1201/ + ^030, 

f2{x,y) = SAsooxy'^ - '2.A2wxy + ^2012/^ + ^1202^ - MiiV + ^021, 
faix.y) = -3^3ooa;^y + A^x^x^ - 2A2oixy + Amx - Aw2y + A012, 
Ji{x) = Azmx^ + A2oix^ + A102X + Aqos- 

The equations in classical mechanics are obtained by setting Ti = 0. We may 
notice from the quantum equations, or directly from the condition [H, X] = 0, 
that we can express all quantum integrable potentials as h^V{x, y) where V 
does not depend on h. It is often more natural and interesting though to choose 
arbitrary parameters contained in V to be depending on ?),, so that the potential 
V does not vanish in the classical limit. One may always verify though that 
through an appropriate transformation of the arbitrary parameters one can write 
the potential as being proportional to 

The three last equations, identical in the quantum and classical cases, yield 
a linear compatibility condition for V , which reads 

= — fsVxxx + (2/2 — 3/4) Vxxy + (~3/i + 2/3) Vxyy " f2Vyyy 

+ 2 {f2y — fzx) Vxx + 2 {~3fiy + f2x + fiy ~ 'ifix) Vxy + 2 {— f2y + fzx) Vyy 
+ (^3/ij,j, + 2/2x1/ — fsxx) Vx + {^f2yy + ^.f^xy ~ 3f4xx) Vy . 

(6) 

Further nonlinear compatibility conditions can be obtained from Q to JSJ 
for the potential, and these are listed in ■ They are quite complicated though 
and were not used for the results stated in this paper. 
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Ill Potentials separable in cartesian coordinates 

If we set V — Vi (x) + V2 (y) in equations 121) to (0) , we find 



^2 

— giVlx + 52V2J, — I flVlxxx + f4V2yyy 

4 V (Y) 
+ 8A3oo(a;V2y - yVi^) + 2 iA2iaVix + ^201^2^) j , 

{9i).=^hiy)Vix + f2{x,y)V2y, (8) 
{92)y^ f3ix,y)V,x+3Uix)V2y, (9) 
(51). + {92). - 2 (/2(x, + Mx,y)V2y) , (10) 

with ?i = in the classical case. Equations © and © are readily integrated, 
so in the cartesian case two equations remain to be solved. 

The compatibility condition © allows us to find ODEs for Vi and V2 ■ If we 
set alternatively y = and x ~ 0, we find 

wx + Aiii)V('{x) + 12yl2ioK(2;) ^ ax + b 

(11) 

(Aaoiy^ - Amy + Ao2i)V^^\y) + 4(2A2oiy - Aiii)V^'{y) + 12A2oiV^{y) ^ cy + d 

(12) 

The solutions to the homogeneous part of these equations are easily found 
and brought to a simple form by translations in x and y. If we take the first 
one for definiteness, we have four different types of solution. When A210 7^ 0, 
we have two possible solutions, 

Vlhon 



^Ihom 9 ^ q 

X X^ 



{x + a)^ {x — a)^ 

Cl . C2 



If yl2io = and Am ^ we get 

Vlhorn = ^ + C2X. 
X'^ 

Finally if only A012 7^ 0, the solution may be brought to the form 

Vlhom = C2X^ + CxX. 

Special solutions are also simple. If A210 7^ 0, we have V\part = ol^^ + P^. 
Otherwise, when Am 7^ 0, Vipart = o:x^ + f3x^, finally, if only Aqu 7^ 0, 
Vipart = ctx"^ + (3x^ . Provided that Hll|l or H12|) do not vanish trivially, we can 
choose Vi or T^, respectively, amongst the following functions: 

(A.l) /i = + + c^x^ + CAX 

(A.2) /2 = ft + fl + C3X2 + C4X 
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Superintegrable potentials 


Leading-order terms 
of the integrals 




^ 1 L ' Jr'XJr'y J 1 L ' Jr^y J ' L ' Jr^X J 


-kVh =a(x^ + y^) + -^7 + 


I 1 Jr'XJr'y j 


Vr =a(x^ +v^) + K + K 


Ij \ Tj Drr.'Di.r 
^ 1 Fxyy J 




1 L ' fxtr'y } 1 L ' J 


V (a;2 + „2N _^ _^ ft^ 






2i3-3a2{L,p2} 


+y ) + F + T^W + 


^9 =^ix +y ) 

. . . . rY 


or'^ o 2/'r7" 21 1 fT" 2l\ 

- 3a^({i,p^} + {L,p;}) 


(y-aY ' {x-aY (y+aV (x+a)^ 


-i-Vi, — of4T^ -1- -1 4- rr 


2 

FxPy 




{L,Pi} 






^fc— gQ4(9a; + y ) + (2^+j^)2 1 (y_Q)2 


{l,pI} 


Vi =^ + 


{L'^,Px}; {L,pxPy}; pI 


Vm — -f 


L3; {L\px}; {L\py} 
{L,pxPy}; pI; pI 




{L,pI}; pI; PxpI 


Vo + V{x) 





Table 1: Superintegrable potentials that satisfy linear compatibility conditions 
for nonzero parameters. 



(A.3) h = ^ + C2X^ + c^x^ + CiX 

(A.4) /4 = CiX^ + C2X^ + CzX 
(A.5) /5 = Cix'^ + C2X 

(A.6) /e = cix2 
(A.7) /7 = cia:, 

and then solve |7J to (|10|l . These long but rather straightforward calculations 
yield the 15 superintegrable potentials included in Table 1. Some of them are 
obviously particular cases of others, but we listed them separately to account 
for their additional integrals. Only the third-order integrals are listed, some of 
them being trivial consequences of lower-order ones. With the exception of the 
harmonic oscillator, potentials that have first-order integrals are not listed here 
for they were already presented in with all their third-order integrals. The 
complete integrals of motion can be found in Appendix I. 

Many of these potentials were not known. The only classical potentials 
among these are indicated with a These are well-known superintegrable 
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potentials (see e.g. [3]), and all of them, except Vi, are in fact quadratically 
superintegrable. 

All the potentials are superintegrable in the quantum case. We therefore no- 
tice that classical nontrivial potentials can have many different quantum equiv- 
alents. The classical harmonic oscillator Va can be seen as a limiting case of the 
quantum potentials Va, Vc and Yd, and also, if we set a = Vft/to, of Vf and 
Vg, not to mention the similar potentials that can be obtained by permutations 
of X and y. The anisotropic harmonic oscillator with ratio 1 : 3 also admits 
many quantum deformations but, interestingly, the anisotropic oscillator with 
ratio 1 : 2 does not admit such deformations. Notice also that if we want to deal 
with real potentials only, a must be either real or purely imaginary in potentials 
Ve, Vf, Vg and Vk ■ Therefore these have as a classical limit harmonic oscillators 
with a > 0. 

All quantum superintegrable potentials reduce to classical ones when the 
classical limit is considered, sometimes in more than one way. For example, 
potentials Vg, Vf and Vg give the free motion potential instead of the harmonic 
oscillator if a remains constant as ?i ^ 0. 

These potentials all satisfy the linear equations and H12|l . and can be 
expressed as sums of simple superintegrable potentials. 

Let us now set ^210 = ^111 = ^012 = so that l(TT|l vanishes trivially. We 
may also assume that Vi does not take one of the forms (AQ to (ACJ, for we 
have already worked these cases out. This is quite useful, for if we set y — 1 
in ©, we obtain for Vi an equation of the same form as with different 
coefhcients. These coefficients must therefore vanish, so A^qo — ^201 = ^102 = 
0. This is a significant simplification that allows us to restrict our attention, 
when considering equation H12|l . to the following three cases: 

i) V2 = ay2. 

ii) V2 = ay; 

iii) A120 = ^021 = 

Before we consider each case separately, it is worth noticing that potentials 
of the form V = Vi (x) that admit third-order integrals independent of y and py 
should appear as solutions here, for the integral remains if we add a function 
^2(2/) to these potentials. These potentials were found in and to satisfy 
equation 

hW;^ = 4Vf - 32^1 - 93, (13) 
and can therefore be written as 

Vi - n^vix), (14) 

where 'P{x) is the Weierstrass elliptic function. Since the y variable plays no role 
here and these potentials admit integrals with leading-order terms proportional 
to , these solutions will appear in cases 10) to (|nHl . 
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A. CaseS V = Vi{x) + ay"^ 

When V — Vi{x) + ay^ and a 7^ 0, wc find that A^qq = 0, and the following two 
equations must be satisfied: 

= Ao3o (?iV/'^ - 6{V,^Y) + 71 y/ (15) 

= A120 (-riV/^^ - 2Aa{xViy + 6{V^)" - Aax^Vl' + 8a^x^^ ^^^^ 

+ 8aAo2i i2ax - (xVl)' - 2Vl) + 4ry(2a - ¥{'). 

where 71 and rj are arbitrary constants. When Ao30 7^ 0, equation p5|l is 
equivalent to (|13|l (up to a translation of Vi to get rid of 71), hence its solutions 
are of the form p4|l . These potentials cannot satisfy simultaneously equation 
H16() for nontrivial parameters. This can be observed by expanding (|13|l in series 
around x = and substituting the result in (|16l) . Therefore solutions given by 
^030 7^ are of no special interest here. 

Let us now set ^oao = 7i = 0. Equation (|16l) can be greatly simplified. 
We assume that 7^ 0, for otherwise equation Hlt)|l can be solved to give 

potential Then by an appropriate translation of x and V , we can get rid 

of the terms involving ^021 and 77 and finally divide by A120; 

= -r?vl^^ ~ 2^a{xVi)' + &{V^)" - Aax^VC + 8aV. (17) 
This equation admits a first integral, namely 

k = h^ {xVi" - y/') + ^x {ax^ - 3Vi) Vi + 6Vf + l2ax^Vi ~ 2a^x^. (18) 

Both (inil and (HHJ can be simplified by setting Vi ^W{x)+ ax^ ji. Then 

^i^vi^W ^ YiWW" + \2{yV'f + hxW' + 2hW - ^6^^ (19) 
with h = -8a ^ for and 

fc2 = {xW" - W") - l8x{W'^)' + 2{2ax^ + iWf. (20) 

for the first integral. Equation p9|l is well known. It is equivalent to equations 
(3.16) in j5j and (2.17) in [20], which were obtained by nonclassical reduction of 
the Boussinesq Equation. It was also shown in [3] to be a nonclassical reduction 
of the Kadomtsev-Petviashvili equation. It has the Painleve property, and, when 
b ^ 0, its solution, given in jS| (equation 2.88) may be written in terms of the 
fourth transcendant function of Painleve, namely 

w = ^b^piix, A) _ ^bP!{x, A) _ hxp.ix, A) _ 1 (^^2 ^ f^2j^^ _ (2^) 
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where bi = ±\/^ = ±\^^8a and P4{x,^) = P^ix, Ki, K2) is the fourth 
transcendant function of Painleve, and therefore satisfies equation 



p;(^, a) - '^P,{x, af - 2axP,ix, af 

-(^a;2 + Xi)P4(x,a) + . 



(22) 



2 ^' -^y ^ ' Pi{x,ay 

Ki and K2 are integration constants. The potential therefore reads 

V{x, y) =a(a;2 + y^) + ^61^^(2;, ^) + AaP^x, ^) 



4aa;P4(a;, — ^) + + hhi). 



(23) 



This potential admits as special cases two anisotropic harmonic oscillators, 
V = aix"^ + 2/2) when K2 ^ (and P4 0), and V = a(x2/9 + y^) when Ki = 
and i4r2 = —1/18 (and P4 = — x/3), as well as all their quantum deformations 
that have the form V — a{p'^x'^ + y"^) + f{x), that is potentials V^, T4, Vj and 
Vfc (up to a permutation of x and y) . 

The constant term in the potential, {—Tt?Ki + fibi), can be set to zero, but 
we will keep it in order to be able to write the quantum and classical integrals 
in a unified way. 

In classical mechanics, the equation (|18|) with h — admits a first integral, 
which reads 

_ (9^1 - ax^){Vi - ax^f + x - ^(^1 - ax^){'iVi + ax^) 

x"^ 

We may therefore write the solution for Vi implicitly as 

cx^ -(f + 2d{Vi - ax^){Wi + ax^) = {9Vi - ax^)iVi - ax^)^, (24) 

where c and d are arbitrary constants. If c = d = 0, we find either the familiar 
anisotropic harmonic oscillators, or a potential obtained by joining at a; = two 
halves of anisotropic harmonic oscillators with different ratios. Even though 
this potential does not have a continuous second derivative, it can be obtained 
as a limiting case of the family of smooth superintegrable potentials 1)26(1 . 

In the general case the potential may be expressed as the root of a fourth- 
order polynomial with three arbitrary parameters (a, c and d), although we can 
set a = 1 and one of the other two coefficients to ±1 by scaling x and Vi. 

Since equation H24() may describe new two-dimensional potentials with bounded 
motion, it is worth studying for interesting special cases. Indeed, if we assume 
a > and d > 0, we can consider the case c = I28a^d'^ /27'^ and d = Aa'^d'^ /27. 
Four solutions exist in that case, two of which are translated harmonic oscilla- 
tors. The remaining two potentials are 



V ^ay^ + ^{2d + 5x^±'ix\/{d + x^)] . (25) 
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Up to an additive constant, potentials (|25|l are equal to 



V = ay^ + ^ (^x±2\/d + x^^ . (26) 

Those potentials are smooth interpolations between anisotropic harmonic 
oscillators with ratio 1 : 1 and 1 : 3. When d — (and therefore c = d = 0) they 
are the junctions of the two halves of harmonic oscillators mentioned above. 

The integral of motion is similar in quantum and classical mechanics, and 
reads 



X = {L,pI} + {ax^y-3yVi{x),p,} - — {— l/j,,, + {ax^ - 3Vi)Vi^, Py}, (27) 

where Vi is a solution to equation l|24|l . Therefore the integral for (|26|l must be 
slightly modified to take into account the constant term we removed. 

B. Caseini V = ay + Vi{x) 
Here we find again two equations: 

= h^Ai2oVl^^ - 6Ai2o{V^^y + 12a2v4oo3 + liV(, 

= ^ (-6(^1^)' + n'vl'^)'- ^Ai2o (6(0:^1)' + x'v;') (28) 

- aAo2i iixVlY + 21//) - 72 + a7i- 

This time we cannot treat the equations separately, but we can, when A120 ^ 
0, translate V to annihilate 71. Then we can substitute the first equation in 
the second to get rid of ^030, and finally translate x to get rid of ^021- We can 
then solve the remaining system, first by solving the second, linear equation, 
and then by substituting the result in the first one. The only solution remaining 
is then potential Vn = ay + Hence we can set A120 = 0, and therefore also 
^003 = 7i = (as Vi 7^ bx and a =/= 0). In the second equation we can set 
72 = 0. If then Aq3o = 0, wc find Vi — ajx^ which was already classified. If 
^030 7^ and Aq2\ = 0, though, the potential Vi{x) is solution to 

?iV/' = 6lf + Xx + k (29) 

If A we can set fc = 0, and we find that Vi can be expressed in terms of 
the first Painleve transcendent: 

Comme A est une constante arbitraire, il est plus pratique d'crire uj^ — 

V = ay + h^uj^Pi (ujx) , (30) 

with = ^ . 

The integral of motion is 
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Notice that in order to consider the hmiting case a = 0, we can multiply the 
integral by a, and thus we find a potential that depends only on x with a trivial 
Py integral. We can also set u)^ = aw^ first, in which case we find the potentials 
l|14(l . If we look for the classical limit, we find 

V = ay + h^/x. (31) 

with the integral 

352 

X = 2pl + 3&{Vi,p4 - {^.Py}- 

Returning to equation (|29|l and assuming A = 0, we find again that Vi has 
the form p4|l . and the integral depends only on x and Px- 

Finally, if we set v4o3o^o2i 7^ 0, we have to solve the equation 

0=(-^(V^i7 + fv^/^^)' + MW + 2y/), 6=^^^0. (32) 
It can be integrated once to give 

Ci - -\{Vl)' + ^V/^^ + h {{xVD + 2V,) . 

We can set Ci = by translations in x and V . The equation admits a first 
integral which reads 

2bn^ {Vi{x) - hx) Vi\x) + hh^ {2b - Vl{x)) Vl{x) - 8bVi{x) {Vi{x) - bxf = ki. 

(33) 

In the classical case, ti = 0, this is enough to solve. The solution, which can 
be written implicitly in a more compact form, is given by 

d ^ Vi{Vi ~ bx)^ , (34) 

where d is an arbitrary constant. When c? = we find the familiar case V = 
bx + ay which admits a first-order integral. Wc may notice that the implicit 
form of the solution is somewhat similar to (|24|) . 

In the quantum case, we notice that the transformation W{x) = ^^1(2;) — bx, 
that preserves the Painleve property, simplifies equation H33(l that becomes 

{2WW" -W'^) -8{W + bx)W^ ^ k2. (35) 

or 

14^" = ^ + !^ + ^ + ^ (36) 
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We can also substitute Y{x) = y^W{x) to find 



If = 0, we can set h = 1 and 6 = 1/2 by the change of variables 



Y = i2hb)^Z, ^=[^J ?■ 

The solution for Z{£^) is then a special case of the second Painleve transcendant, 
defined by the equation 

P^'ix, a) = 2P2ix, + xP2{x, a) + a. 
The solution for V is 

V{x,y)=bx + ay+{2hb)ip^ (^l^^y x,Q^ , (37) 

with the particular case V = bx + ay. 

Now if fc2 7^ 0, we use . that can be normalized by a change of variables, 



We then find 



Y" ^^+APY'-^Y-^, (38) 

where /3 = — \/— fe/(4?i&) 7^ is an arbitrary constant. Equation H38|l corre- 
sponds to case XXXIV, p. 340, in ^^l- The solution for F(^) reads 

2(3Y = P^iC, -2(3 - i) + P2(^, -2/3 -\f + \ (39) 

where P2 is once again the second Painleve transcendant. Since /3 is an arbitrary 
constant, we can set n = —2(3 — i. 

Back to the original variables, we get 

V{x,y)^ay+{2n'b'Y (pl^i-iAb/fi^)"^ x, n) + Pl{-{Ab/h^f^ x, k)) (40) 

This potential admits V — ay and V — ay + T? jx^ as particular cases. 

In all cases with V — ay^ ^1(2;) ^J^d ^012^030 7^ the integral of motion is 

2avl - 2bplpy + a{Wi{x) - bx,p^} - 2b{Vi{x) , py} . (41) 

Limiting values for a, b and h give either known potentials or trivial integrals. 
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C. CaseED ^120 = ^21 = 



Since we set here ^120 — ^021 — 0, we assume ^030^003 7^ (otherwise we 
would find potentials of the form 114|) '). The conditions for the existence of a 
third-order operator then read 

h^V{'{x) = 6V^{x)+Aoo3CTx 

(42) 

h^V^'{y)=6Vi{y)-Ao3o<yy- 
If a vanishes, we find the potential 

with two integrals that each depend on only one variable, as could have been 
predicted from the results of equation (|13|l . If a does not vanish, let us set 
61 = ^0030' and &2 = -Aqsoct. 
If ?i = 0, we find that 

V = ±^/(3^ ± ^/(hy 

is superintegrable, where the Pi = bi/6 are arbitrary constants. 

This defines similar real potentials on each quadrant of the plane. We can 
again patchwork the pieces to find real continuous potentials defined everywhere, 
e.g. V = ci \/\x\ + C2 \/[yf, although in that case neither the Hamiltonian nor 
the integral are differentiable at the equilibrium point. 

Finally, if ?i 7^ 0, we find that Vi and V2 can be both written using the first 
Painleve transcendent. 

V = h^LulPi {toix) + T^ojIPi {uj2x) (43) 

where the uoi — (&i/?i^)^/^ are arbitrary constants. 
The integral of motion is 

X = 262pi^ - 26ip3 + 3&2 {14 (a; , 61 ) , p J - 36i { F2 (2/, 62 ) , P J , (44) 
both in quantum and classical mechanics. 

IV Conclusion 

We have found all systems in two-dimensional euclidian space that admit sepa- 
ration of variables in cartesian coordinates and at least one third-order integral, 
both in quantum and classical mechanics. Many new superintegrable potentials 
were found, and, interestingly, all the quantum superintegrable potential are 
found as solutions of equations having the Painleve property, and this is prob- 
ably not accidental. Many of the quantum potentials can in fact be written in 
terms of different transcendent functions of Painleve, and many are related to 
group invariant solutions of the Korteweg-de Vries equation, and to reductions 
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of the Boussinesq and KP equations. All classical integrable potentials were 
found to be limiting cases of quantum ones. They do not obey to equations 
having the Painleve property, though, for many classical integrable potentials 
have movable branch points of the form \Jx — h. Thus, in that respect, quantum 
integrable potentials behave more regularly than classical ones. A natural ques- 
tion is what does the Painleve property, in quantum mechanics, tell us about 
classical integrable potentials. Since this paper deals mostly with the classifica- 
tion of superintegrable systems, the consideration of such questions regarding 
their properties and their solutions will be postponed to a future article. 

Our investigation provided interesting new examples of the differences be- 
tween quantum and classical integrability. A systematic search for systems with 
higher-order integrals is therefore a useful task since we still know little about 
such systems, and they are likely to share interesting properties. 

Note added in proof: The future article mentioned earlier (55) indeed shows 
that superintegrability of separable systems tells a lot about the physical prop- 
erties of the separated systems. Even though the systems considered here can 
all be separated in two independent one-dimensional systems, two-dimensional 
superintegrability provides information on them that is far from obvious from 
a one-dimensional perspective. 

Equations (01 to lO are the necessary and sufficient conditions for the exis- 
tence of third-order integrals. The general solution to these equation is highly 
nontrivial, and it is not likely that a direct approach will lead to such a solution 
for integrals of order higher than four, without the use of new methods. In 
order to develop such methods, it would be useful to have a rigorous definition 
of quantum integrability. 
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V Annex I 

Here is the complete list of two-dimensional Hamiltonians separable in cartesian 
coordinates that admit at least one third-order integral. We also gave all their 
third-order integrals. Potentials depending on only one cartesian coordinates 
are not listed here, for they were classified in 12 . Many of the potentials 
listed below were already known to be superintegrable, but we listed them here 
for completeness and in some cases to take into account additional third-order 
integrals. Potentials in boxes are those who do not admit enough first- or 
second-order integrals to make them superintegrable. Most of them were not 
known before. 

Some potentials (namely and Vq) are not superintegrable according to 
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the usual definition, for their integrals have simple though nontrivial relations 
with the Hamiltonian. For definitions of parameters the reader is referred to 
the body of the article. In order to write the integrals of motion in a compact 
form we often use the notation V{x, y) = Vi{x) + V2{y) = Vi + V2. 



A. Quantum potentials 

(Q.l) y = a(a;2+y2) 

X2 = {L,PxPy} + a{2x'^y,py} - a{2xy'^,Px} 
X3 = {L,pI} + 2a {{xy^,Py} - 

= {L,pI} - 2a {{x'^y,Px} - {x^,Py}) 

(Q.2) V = a{x^ + y^) + X + ^ 



X-L = {L,pxpy} + {xy (-|| + 2ax) ,py} - {xy (^-p + 2ay^ ,Px} 



(Q.3) y = a(x2+ 2/2) + 1^ + 1^ 

X, = 2L3 _ ?,2|3^ + 2y + + +2x+ ^-f,py} 

X2 = {L,pa,py} + {xy (-^ + 2ax^ ,py} - {xy (-f^ + 2ay^ ,p^} 

(Q.4) V = a{x'+y^) + ^ 

X,=2L^-n\2y+^-^,Px} + n\^4- 



2X,Py} 
{ 

X3 = {L,pI} - {2ay^ + n^,p,} + {H^^ + 2axy^py} 



X2 = {L,PxPy} - 2{axy^ - \f,Px} + 2a{x^y,py} 



(Q.5) 



v=n- 



[x — a)"^ [x + a)2 



X, = 2L^ - ia\L,pl} + i^{yi-8 + ^ - (i^fet^),^.} 



4 a2{x-a)2{x+a)2 JiPy/ 



4a4 



X2 = {L,pl} + h'{yi^'^ - 



-2 J- 2(2 -4a ) 

27? 



2x , i.r(.i-+<x-) 

(x-cx)-'{x+cxy^'i'yi 



(Q.6) 



1 



1 



(3; + q:)2 (a; — a)^ 



(x— a)2(a:+a)^ 
(,r-n)2(.r + ,>)2 >■ I'tiS 



(Q.7) 



1 



+ 



(y — a)2 (x — a)2 {y + a)^ (x + a)^ 



X, = 2L' - 3a^{L,pl} + {L,pl}) + f {y(124 + ^ii-jU + 



2'l{x^-5y'') 
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a;2_Q2 -t- (a;-Q)2(2; + Q)2 ^ (y-Q)2(y + Q)2 ^ -Pk J" 4\-^V-L^^"r „2 

24(5a:''-y'') _ 144a:'' _ 24(x^-3j/^)(j/^+a^) , 48(a''-j/'')(s'' W) ^ 



(x— a)^(x+a)^ 



(Q.8) 

(Q.9) 

(Q.IO) 
(Q.ll) 



Xi = 2p^pl + {-2ay2 + p^^^} + {gaxy + cy,Py} 



Xi = {L,pI} + ^a{y^,p^} - 6a{xy^,py} 



V = ai9x'+y') + ^ 



X, = {L,pI} + - f,pA + {^x (-2ay2 + ^) ,py} 



V=h'(^{9x'+y') 



1 



+ 



X,^{L,pl} + n'{y 



8a^ 



(Q.12) 





a 


x-^ 





Xi = {L\p,} - 2h^{^,py} + {3h% + 2a^ 
X2 = {L,p^py} - 2h^{^,py} + 2a{^,px} 
X3 = 2pl + {^,p,} 



(Q.13) 



V = 



X2 = {L\p,} - 2n^{lpy} + {3n% + 2T^^ 



(Q.14) 



Xs = {L\py} - 2n\^,p,} + {2.%^^ + 2n% 

Xi = {L,P,Py} - 2h\^,Py} + 2h\^,P,} 

X, = 2pl + {^,p^} 
X^ = 2pl + {^-f,Py} 

V = ax+^ 

X, = {L,pI} - {f + {3^ - a^,p,} 

X2 = 2p3 + {3|^,pJ 

^3 = 2p^pl + 2te{jj,p^} + a{y,py} 
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(Q.15) V = n''V{y) + V{x) 



Xi = 2pl + {3n'V{y),Py} 



(Q.16) V=h\r{x)+P{y)) 



Xi = 2pl + {3h^V{x),p,} 
X2 = 2pl + {3h^V{y),Py} 



(Q.17) V= n^JlPi {ujix) + n^ujlPi {uj2x) 



2, ,2 1 



Xi = 2ojIpI - 2uj\pI + 2,uj^{Vi{x,uji),p^} - 'i^^\{V-2{y,uJi),Py) 



(Q.18) 



V = a{x^ + y^) + ^biPi{x, + 4aP^{x, + AaxP^ix, + ^{-Ti^Ki + hbi) 



Xi = {L,pI} + {ax^y - SyVuPa:} - + {ax'' - SVi)Vi^,Py} 



(Q.19) V = ay + h'co'Pi{ujx) 



Xi = 2pl + 3nV{Pi(u;x),p,} + {^,Py} 



(Q.20) 




Xi = 2apl - 2bplpy + a{3Vi{x) - bx,p4 - 2b{Vi{x),py} 



(Q.21) 



V = ay+ {2n%'') ^ (^Pi{-{Ab/n'')ix, k) + P^{-{Ab/%')^x, k)) 



Xi = 2apl - 2bplpy + a{Wi{x) - bx,p^} - 2b{Vi{x),Py} 



B. Classical potentials 

(C.l) y = a(.T2 + y2) 
Xi = L3 

X2 = {L,p^py} + a{2x''y,py} - a{2xy'',p^} 

X3 = {L,pI} + 2a {{xy',Py} ~ 

X4 = {L,pI} - 2a {{x'y,Px} - {x^,Py}) 

(C.2) V = a{x' + y') + ^ + ^ 



Xi = {L,p^py} + {xy (-|| + 2ax) ,Py} - {xy (-p + 2ay) ,p^} 

(C.3) V = 0(4x2 + y-^) + ^+cx 

Xi = 2p^pl + {-2at/2 + + {8axy + cy,Py} 



(C.4) y = a(9a;2+y2) 



Xi = {L,pI} + |a{y^p^} - 6a{xy^,py} 
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(C.5) y=±VA^±V^ 



(C.6) ' ' ' 



V=ai/ + Vi 


where Vi satisfies equation (|24|l 
{ax'^y - 3yVi,Pa:} - j^{{ax'^ - iVi)Vix, Py} 


x,^{l,pI} + 


V= ay + h\fx 




V= ay + Vi{x 


) , where (Vi - bx)^Vi ^ d 



Xi = 2apl - 2bplpy + a{3Vi{x) ~ bx,p^} - 2b{Vi{x),Py} 
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